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STOCHASTIC  MODELS  OF  THE  DISTRIBUTION  OF  DYADIC  WARFARE  IN  TIME 


Stochastic  models  are  constructed  to  illuminate  the 
dynamic  incidence  of  international  warfare  during  the  l8l6- 
1965  period.  It  is  argued  that  the  probabilistic  structure  is 
revealed  most  clearly  through  an  analysis  based  on  dyads  of 
nations,  thereby  disassembling  such  large  wars  as  World  War  II. 
The  conceptual  focus  is  maintained  on  a clear  delineation  of 
heterogeneity  over  time  and  over  actors  and  of  contagion  in 
both  its  addictive  and  infectious  varieties.  Departures  from 
randomness  are  considered  as  modifications  of  the  Poisson  process. 
Methodological  attention  is  directed  at  statistical  analysis  of 
the  interarrival  times  between  initiations  of  dyadic  warfare. 
Cyclic  behavior  is  investigated  tnrough  a cosine  wave-form  vari- 
ant of  the  Poisson  process.  A conclusion  of  infectious  behavior 
is  supported  by  a variety  of  analyses.  An  autoregressive  model 
of  order  4 is  found  to  adequately  fit  the  interarrival  times  and 
account  for  the  infectious  behavior. 


1.  THE  PROBLEM 


The  initiation  of  warfare  is  a dramatic  event  on  the  world 
scene,  absorbing  attention,  resources,  and  lives,  often  on  a 
grand  scale.  War  , in  the  conventional  Western  view,  is  a cata- 
clysmic breakdown  of  the  international  order,  a catastrophe; 
i.e.  war  between  two  nations  becomes  the  ultimate  dyadic  inter- 
action. Its  onset  is  fraught  with  uncertainty  as  the  basic  unpre- 
dictability of  timing,  participants,  number  of  casualties  and 
level  of  destruction  becomes  strikingly  clear.  The  study  of  war, 
in  all  Its  parts,  has  proved  fascinating  for  centuries  and  has 
been  pursued  with  diligence  from  many  different  perspectives, 
including  the  most  formal  and  theoretical. 

Perhaps  the  most  effective  beginning  on  understanding 
the  conditions  of  war  through  the  use  of  formal  mathematical 
models,  tempered  by  statistical  verification,  was  made  by  Lewis 
Fry  Richardson  in  his  path-breaking  work.  The  Statistics  of 
Deadly  Quarrels  (i960)1.  It  has  become  increasingly  evident 
since  the  appearance  of  his  work  that  a particularly  effective 
way  of  coping,  in  a theoretical  sense,  with  the  complex  inter- 
actions and  uncertainties  attendant  to  warfare  is  through  the  use 

of  explicitly  stochastic  models,  formally  treating  uncertainty 

2 

through  probability  theory. 

We  demonstrate  in  this  article  that  some  of  the  fundamental 
structural  questions  that  have  commanded  the  attention  of  theoret- 
ical analysts  can  be  addressed  using  stochastic  models  coupled 
with  statistical  inference  procedures.  In  particular  we  shall 
focus  our  efforts  on  a clear  delineation  of  the  concepts  of 
heterogeneity  and  contagion,  as  useful  theoretical  constructs 


in  a multi-actor,  interactive  system. 
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Heterogeneity  and  contagion  will  be  considered  as  devia- 
tions from  the  randomness  implied  by  the  Poisson  process  model. 
Richardson  (i960,  pp.  129-142)  reported  a good  fit  of  the  distri- 
bution of  wars  in  time  to  the  Poisson  process  model,  using  data 
derived  from  Wright's  (1942)  list  of  wars.  Others,  including  Moyal 
(1949),  Denton  (1966),  Denton  and  Phillips  (1967),  and  Singer  ana 
Small  (1972),  have  also  reported  research  in  this  direction.  Their 
research  has  focused  on  the  alternative  of  a time-heterogeneous 
departure  from  randomness  required  by  cyclic  behavior  in  the 
periodicity  of  warfare.  Richardson  (i960,  p.  173  and  p.  265) 
has  also  evaluated  "war-proneness"  as  one  conceptualization  of 
heterogeneity. 

A conceptually  different  tack  was  followed  by  Singer  and 
Small  (1974)  in  evaluating  the  proposition  that  "war  begets  war", 
a basic  assertion  of  contagion.  While  they  term  this  proposition 
"part  of  the  folklore  of  international  politics",  it  is  nonethe- 
less important,  since  if  contagion  were  a determinative  factpr, 
it  would  offer  a conceptually  important  explanation  of  the  occur- 
rence of  warfare.  Recognition  of  the  influence  of  contagion 
would  then  be  useful  for  prediction  purposes  about  the  growth  of 
warfare  and  motivate  the  examination  of  various  causal  factors 
stimulating  or  inhibiting  contagion  in  warfare. 

To  investigate  the  contagion  proposition  Singer  and  Small 
use  their  data  set  on  wars  between  1816  and  1965  (Singer  and  Small, 
1972).  Using  these  data  allows  them  "to  ask  whether  those  years 
during  which  international  war  began  were  more  likely  to  be  fol- 
lowed by  another  war  in  the  same  or  subsequent  year  than  were  those 
in  which  no  war  began"  (Singer  and  Small,  1974,  p.  279)-  Their 
"war  vs  no-war"  2x2  turnover  table,  however,  shows  a pattern  which 


is  not  a statistically  significant  departure  from  an  independence 
hypothesis.^  Singer  and  Small  thus  conclude  that  war  is  not 
contagious. 

This  is  not  a surprising  finding.  Indeed,  it  could  have 
been  predicted  on  the  basis  of  the  analysis  presented  by  Singer 
and  Small  (1972)  in  Chapter  9 of  The  Wages  of  War,  "Cycles  and 
Periodicity  in  the  Incidence  of  War."  There  it  was  shown  that  the 
empirical  distribution  of  the  intervals  between  the  onset  of  wars 
offered  a good  fit  to  an  exponential  distribution.  The  exponential 
distribution  is  derived  from  assumptions  implying  randomness  and  is 
consistent  with  a Poisson  process  model  for  the  occurrence  times  of 
wars. 

Although  one  basic  conclusion  of  Richardson  (i960)  was  that 
the  onset  of  wars  fit  a Poisson  distribution,  some  of  Richardson's 
analysis  appears  to  contradict  randomness  and  point  in  the  direction 
of  contagion.  This  emerged  out  of  Richardson's  attempt  to  construct 
a mathematical  model  which  would  adequately  describe  the  number  of 
nations  on  each  side  in  a war.  In  constructing  this  model  Richardson 
went  through  a lengthy  process  of  making  assumptions,  forming  these 
into  a model,  and  then  either  discarding  the  model  (for  a variety 
of  reasons)  or  testing  it  against  his  data.  Richardson  considered 

some  thirteen  theories  with  the  twelfth  giving  the  best  fit  to  the  data. 
The  assumptions  which  made  up  Richardson's  theory  XII  were  essentially 
that  the  number  of  parties  on  each  side  in  a war  was  the  outcome  of 
a process  heavily  influenced  by  geography  and  modified  by  infection. 

This  would  appear  to  contradict  Singer  and  Small  as  well  as  Richardson' s 
own  conclusion  of  randomness. 
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The  contradiction,  however,  is  more  apparent  than  real. 

Richardson' s finding  of  infection  was  based  upon  an  analysis  which 
used  the  Individual  nations  as  units,  whereas  the  data  used  by  both 
Singer  and  Small  and  Richardson  which  fit  the  random  pattern  used 
wars  as  the  units  of  analysis.  The  conclusion  which  follows  from 
this  is  that  one  war  does  not  generate  other  wars  (Singer  and  Small, 

1974,  p.  279  > 1972,  chapter  9 ; Richardson,  i960,  pp.  128-31  ), 

but  that  once  a war  starts  the  chances  of  other  nations  being  drawn 
into  that  war  may  be  influenced  by  infection  (Richardson,  i960, 
pp.  275-86). 

It  is  possible,  moreover,  that  Richardson  understated  the 
extent  to  which  wars  were  enlarged  by  the  contagion  of  the  conflict. 
3ecause  his  models  were  exceedingly  complex  and,  by  contemporary 
standards,  his  computational  ai-as  rather  primitive,  Richardson  did 
not  include  in  his  final  theories  (including  XII),  those  wars  in 
which  there  was  a total  of  more  than  four  participants.  This  effectively 
excluded  17  of  the  91  wars  in  his  data  base,  or  about  13.7  per  cent. 

Since  the  17  excluded  wars  are  precisely  those  which  would  have 
exhibited  the  greatest  amount  of  contagion,  it  is  perhaps  now  under- 
standable why  Richardson  found  the  process  by  which  parties  are  added 
to  a war  only  "modified  by  infectiousness." 

The  possibility  of  wars  growing  by  a contagious  process  is 
potentially  interesting,  since  it  could  account  for  the  peaks  of 
war  concentration  which  have  been  noted  by  several  investigators 
(Singer  and  Small,  1972;  Richardson,  i960;  Moyal,  1949). 
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This  paper  reports  an  investigation  of  the  extent  to  which 
wars  are  heterogeneous  or  contagious  in  their  stochastic  develop- 
ment. The  distinctions  between  what  is  reported  here  and  previous 
research  are  several  and  are  important:  First,  the  models  used 

here  are  substantially  different  from  those  used  by  Richardson, 
in  which  no  time  dependency  was  present.  Our  models  are,  as  will 
be  shown,  time-dependent  and  hence  dynamic.  Second,  the  models 
we  will  present  below  are  simpler  than  those  used  by  Richardson. 

Not  only  does  this  simplicity  provide  greater  parsimony  as  a virtue 
in  itself,  but  it  allows  us  to  include  data  not  considered  by 
Richardson.  Third,  the  data  we  shall  use  in  exploring  the  models 
will  consist  of  dyads  which  actually  fought  each  other  in  inter- 
national wars  between  1816  and  1965-  This  is  a longer  time  period 

e; 

than  that  covered  in  Richardson's  investigation.  Fourth,  this 
paper  will  focus  on  heterogeneity  and  contagion  as  conceptual 
models  for  the  initiation  of  dyadic  war. 

* 

In  addition  this  paper  will  explore  and  demonstrate  the 
utility  of  some  methodological  techniques  which  have  not  been 
widely  used  in  political  science.  In  particular,  use  will  be 
made  of  a stochastic  model  of  interarrival  times.  This  is  not 
the  first  time  that  inter-arrival  times  have  been  statistically 
analyzed  in  the  study  of  war  (Singer  and  Small,  1972,  p.  205), 
but  this  is  the  first  time  the  interarrival  times  have  been  ex- 
plicitly modelled. 

We  begin  with  a conceptual  discussion  of  the  nature  of 
heterogeneity  and  contagion  in  a multi-actor  dynamic  system. 
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2.  HETEROGENEITY 

The  processes  of  international  relations,  specifically 
during  crisis  situations,  are  seldom  "well-behaved"--events  are 
complex,  influenced  by  a multitude  of  factors,  and  subject  to  con- 
siderable variability.  This  suggests  that  international  crisis 
behavior,  particularly  during  the  confusion  and  turmoil  of  inter- 
national war,  is  characterized  by  heterogeneity,  rather  than  homo- 
geneity. 

This  heterogeneity,  however,  can  evidence  itself  in  two  quite  dis 

tinct  ways.  First,  the  stochastic  process  governing  the  development  of 
the  crisis  may  be  heterogeneous  over  time,  i.e.,  the  evolution  of  the 
process  may  be  governed  by  different  rules  at  different  times.  This 
time  heterogeneity  may  appear  continuous  in  form,  as  when  the  rate  of 
initiation  of  conflict  accelerates  over  time,  or  abruptly  in  form, 
as  when  the  rate  of  initiation  of  conflict  suddenly  jumps,  due  to,  say, 
an  unexpected  change  in  weapons  technology.  Effects  of  this  nature 
may  be  accounted  for  in  stochastic  models  by  allowing  a process  rate 
parameter,  X,  to  vary  with  time,  t,  perhaps  X = a + bt  o£  by  building 
in  discrete  shifts  in  the  process  rate  parameter,  such  that  the  rate 
equals  \±  and  X2  for  t < tQ  and  t > tQ,  respectively.  There 
is  also  the  possibility  of  cyclic  behavior  of  the  process  rate  parameter 
In  fact,  Denton  (1966)  and  Denton  and  Phillips  (1967)  suggest  a 
cyclical  pattern  of  war  with  peaks  successively  increasing  in  size, 
while  Singer  and  Small  (1972)  find  evidence  to  support  the  conclu- 
sion of  periodic  20  year  peaks  in  warfare.  A basic  question  addressed 
in  this  paper  is  whether  such  observed  behavior  should  be  attributed 
to  cyclical  changes  in  the  process  rate  parameter  or,  rather,  to 
contagious  behavior  as  it  is  delineated  in  Section  3. 
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Second,  the  stochastic  process  may  be  heterogeneous  over 
actors,  l.e.,  the  probabilistic  behavior  of  different  actors  in  a con- 
flict environment  may  be  different  and  not  governed  by  the  same 
probability  laws. 

The  first  type  of  heterogeneity  requires  that  an  adequate 
model  be  time-dependent  in  the  parameters  governing  the  stochastic 
evolution  of  the  system.  Thus  the  parameters  of  the  stochastic 
process  are  made  explicitly  functions  of  time.  As  Ginsberg  (1971, 

1972a)  has  indicated,  some  sources  of  time  heterogeneity,  which  do 
not  interact  with  the  states  or  the  stochastic  development  of  the 
process,  such  as  broad  socio-historical  changes,  can  be  eliminated 
using  the  device  of  operational  time  which  essentially  alters  the 
time  scale  in  a nonlinear  fashion  to  obtain  time-independence  of 
the  parameters.  This  approach  will  be  illustrated  in  Section  5.2.^ 

The  second  type  of  heterogeneity  follows  from  the 
realization  that  we  are  dealing  with  a multivariate  time  series, 
tracing  a cohort  of  actors  as  they  move  stochastically  through 
various  possible  states.  In  international  relations  the  actors  will 
typically  be  nations  or  groups  of  nations  and  hence  will  differ 
substantially  according  to  almost  any  conceptually  valid  set  of 
factors  influencing  the  phenomena  under  study.  Thus  it  may  well  be  an 
inadequate  approximation  to  assume  that  each  of  the  actcrs  is  governed  by 
the  same  probability  laws,  i.e.,  that  the  parameter  values  for  the 
various  actors  are  the  same. 
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CONTAGION:  ADDICTION  VERSUS  INFECTION 

Contagion  refers  generally  to  a variety  of  stochastic 
dependencies  concerning  one  or  more  actors  in  an  interactive 
system.  Unfortunately,  the  term  has  been  used  in  the  literature 
of  stochastic  processes  in  a highly  ambiguous  manner,  so  ambig- 
uous, in  fact,  that  it  has  precluded  adequate  conceptualization 
of  one  of  the  most  important  aspects  of  the  study  of  multi-actors 
systems  (including,  of  course,  those  in  international  relations). 

For  our  purposes  it  is  useful  to  distinguish  between  two 
varieties  of  contagion,  which  shall  be  termed  here,  addiction 
and  infection.  Both  are  relevant  to  the  analysis  of  international 
crises  and  the  growth  of,  warfare.  Addiction  is  a characteristic 
of  individual  actors  and  is  exhibited  in  its  positive  form  when 
the  fact  that  an  actor  has  taken  an  action  makes  it  more  likely 
that  the  actor  will  take  similar  actions  again  in  the  future. 

Thus,  in  international  crises,  addiction  is  evident  when  an  initial 
hostile  act  by  a state  means  it  is  more  likely  that  this  state 
will  accelerate  the  frequency  and/or  hostility  of  subsequent  acts. 
Addiction  may  be  conditional  as  when  aggressors  who  meet  with 
success  become  more  likely  to  be  aggressors  again.  Negative 
addiction  is  clearly  also  possible  since  an  experience  may  have 
an  inhibitory  effect  on  future  occurrences.  Infection,  on  the 
other  hand,  is  characteristic  of  a group  process,  and  is  exhibited 
when  the  fact  that  one  actor  has  taken  an  action  changes  the 
probability  of  a second  actor  taking  an  action  (Coleman,  1964, 
p.  299).  Infection  may  therefore  have  positive  impact  on  the 
behavior  probabilities  of  other  actors  or  negative  impact.  In 
international  conflict  infectious  contagion  is  not  at  all  unlikely. 
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For  example,  the  engaging  in  hostile  interchange  within  one  dyad 
may  significantly  increase  the  probability  that  a third  party 
initiates  warfare  with  one  (or  both)  of  the  nations  in  the  dyad.' 

The  distinction  between  addiction  and  infection  is  important 
on  more  than  conceptual  grounds  for,  as  Taibleson  (1974)  has  noted, 
there  is  considerable  difficulty  in  statistically  sorting  out  whether 
a process  is  heterogeneous  among  the  actors  or  is  addictive.  Taibleson 
suggests  that  to  verify  an  addictive  effect  by  looking  at  a sample  and 
then  looking  to  its  future  to  see  if  prior  events  correlate  with  later 
events,  does  not  determine  whether  the  effect  is  true  addiction 
(where  the  fact  of  an  earlier  event's  occurrence  changes  the  probability 
of  the  occurrence  of  a later  event)  or  spurious  addiction  (where  the 
fact  that  an  event  did  or  did  not  occur  at  an  earlier  time  to  an  actor 
changes  the  estimate  of  the  probability  that  the  individual  came  from 
a lower  or  a higher  risk  stratum  of  the  population).  In  either  case 
you  will  always  get  the  appearance  of  addiction.  This  statistical 
difficulty  is  not  operative  in  a model  assuming  the  infectious  variety 
of  contagion,  provided  data  are  available  on  the  behavior  over  time 
of  the  various  actors.  Since  infection  is  likely  to  be  an  important 
mechanism  of  contagion  in  international  relations  processes,  it  may 
be  possible  to  distinguish  heterogeneity  among  actors  from  contagion 


on  a statistical  basis. 
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4.  DATA 

The  data  set  used  in  this  research  results  from  an  examination 
of  international  warfare  during  the  period  1816  to  1965  inclusive. 

The  basic  unit  of  analysis  is  the  dyad, composed  of  two  nations  which 
have  entered  into  armed  conflict.  The  use  of  dyads  appears  to  be  an 
appropriate  device  for  our  interests.  There  are  several  reasons  for 
this:  Most  importantly,  recording  the  day  on  which  a particular 

dyad  began  fighting  provides  the  basic  time  record  for  evaluating  the 
stochastic  models  incorporating  heterogeneity  or  contagion.  Specifically 
it  allows  us  to  explore  heterogeneity  over  time  and  actors  as  well  as 
the  contagious  concepts  of  addiction  and  infection.  Further,  since 
the  number  of  dyads  in  a war  is  one  surrogate  for  its  size,  it  will 
reflect  the  growth  over  time  of  those  conflicts  which  do  enlarge.  From 
this  perspective  World  War  II  was  not  a single  event  which  started  in 
1939,  but  a gradual  process  which  grew  from  1939  to  1945.  Also  looking 
at  dyads  provides  a more  realistic  approach  to  large  wars,  which  in 
important  respects  were  a complex  collection  of  deadly  quarrels  rather 
than  the  confrontation  of  monolithic  coalitions  depicted  in  some  accounts 

Our  focus  upon  dyads  is  also  based  upon  a definition  of  war  as 
Joint  belligerent  activity,  that  is,  both  nations  must  fight.  The 
wisdom  of  this  definition  may  not  be  self-evident:  There  may  be  a 

tendency  to  think  that  if  one  nation  attacks  another  that  defense  or 
counter-attack  is  automatic,  therefore  suggesting  an  analysis  based 
on  the  individual  nation.  We  would  only  agree  that  there  is  a high 
probability  of  defense  or  counterattack,  while  pointing  out  that  history 
contains  numerous  examples  of  one  nation  launching  an  attack  and  the 
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other  nation  involved  not  resisting.®  No  war,  as  we  typically  view 
it,  existed  in  these  cases. 

To  create  our  data  base  we  have  drawn  upon  the  listing  of 
international  wars  between  1816  and  1965  compiled  by  J.  David 
Singer  and  Melvin  Small  in  their  encyclopedic  The  Wages  of  War  (1972). ^ 
The  basic  advantage  of  using  this  list,  as  mentioned  above,  is  the 
careful  and  rigorous  work  of  the  investigators  in  compiling  the  data. 
Notwithstanding  th'  availability  of  the  data  and  the  care  exercised 
in  its  production,  it  was  necessary  to  make  some  adjustments  in  the 
data  to  make  it  suitable  for  our  purposes.  These  changes  were  largely 
occasioned  by  a central  interest  in  the  dyads  which  actually  fought  each 
other  in  the  war;  we  are  not  interested  in  those  dyads  in  which  there 
is  no  declaration  of  war  or  no  actual  combat.10  This  decision  neces- 
sitated a scrutiny  of  Singer  and  Small's  list  of  opponents  to  determine 
who  actually  fought  whom.  Alterations  were  necessary  only  in  the  cases 
of  three  wars:  The  Seven  Weeks  War  (1866),  World  War  I and  World  War  II. 
In  each  case  the  number  of  dyads  was  reduced  with  the  largest  reduction 
in  World  War  II.11 

The  day  on  which  each  of  the  dyadic  war  initiations  took 
place  in  the  150  year  period  was  recorded.  If  these  data  are 
aggregated  over  the  period  of  a year  (as  in  Table  1),  there  were 
95  years  with  no  dyads  initiating  fighting,  27  years  in  which  one 
dyad  initiated  fighting,  8 years  with  2,  6 years  with  3,  3 years 
with  4,  2 years  with  5,  10,  and  14,  and  1 year  with  6,  7,  11,  13, 
and  40.  This  produced  a total  of  209  dyadic  war  initiations. 

Much  of  the  analysis  in  the  following  sections  will  be 
based  on  the  times  between  successive  war  initiations  (the  inter- 
arrival times).  The  time  between  the  ith  and  (i+l)st  war  initiation 
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will  be  denoted  by  T^.  A problem  which  exists  in  tne  data  is  that 
the  exact  interarrival  times  are  not  available  since  the  times 
are  not  recorded  in  hours,  minutes,  etc.  In  any  case,  recording 
the  data  more  accurately  than  the  day  of  war  initiation  would 
lead  to  essentially  the  same  conclusions.  In  recent  dyadic 
war  initiations,  it  would  be  possible  to  ascertain  the  exact 
moment  hostilities  commenced.  However,  this  is  the  exception 
as  exact  times  could  not  be  found  for  early  wars. 
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5.  DATA  ANALYSIS 

P 

This  section  uses  the  interarrival  times  and  the  aggregated 

data  of  Section  4 to  discuss  the  heterogeneity,  contagion,  and 

randomness  of  dyadic  war  initiations.  Since  the  hypothesis  of  a 

random  distribution  of  warfare  in  time  has  been  supported  by 

Richardson  (i960)  as  well  as  by  Singer  and  Small  (1972  and  1974),  it 

is  appropriate  to  begin  by  assessing  the  extent  to  which  these  data 

fit  the  Poisson  distribution,  a distribution  consistent  with  randomness 

12 

and  inconsistent  therefore  with  contagion  or  heterogeneity. 

5-1  The  Poisson  Process 

.<*'■  Three  assumptions  are  basic  to  and  imply  the  Poisson  process: 

(1)  The  number  of  dyadic  war  initiations  in  one  time  period  is  a 
random  variable  independent  of  the  random  number  of  dyads  initiating 
war  in  another  (non-overlapping)  time  period. 

(2)  In  a sufficiently  short  time  period,  the  probability  of  two  or  more 
instances  of  dyadic  war  initiation  is  negligible. 

(3)  For  sufficiently  short  time  periods,  the  probability  that  a dyad 
initiates  war  is  proportional  to  the  length  of  the  time  period. 

These  assumptions  lead  to  the  following  mathematical  model 

# 

for  the  probability  of  k(k=0, 1, 2, . . . ) instances  of  dyadic  war 
initiation  in  the  time  period  from  tQ  to  t^: 

e-X(ti-t0)  [X(tl -t0)]kAl, 

where  e is  the  base  of  the  natural  logarithm  and  the  parameter  X.  > 0 
is  the  mean  number  of  instances  of  dyadic  war  initiation  in  a time 
period  of  unit  length  (X  is  also  called  the  rate  in  the  sequel). 
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One  of  the  characterizing  features  of  the  Poisson  process 
is  that  the  interarrival  times  between  successive  occurrences  are 
independent  and  identically  distributed  random  variables  with  an 
exponential  probability  distribution.  Thus  if  we  let  T be  the 
random  variable  giving  the  time  interval  between  one  dyadic  war 
initiation  and  the  next,  T would  have  an  exponential  distribution 
with  mean  \ . The  probability  density  function  of  T would 

then  be  given  by  fT(t)  = \e"^,  t>0.  Some  insight  to  the 

implications  of  this  fact  about  the  Poisson  process  is  provided  by 
the  following  derivation:  View  T as  the  random  variable  giving 

the  waiting  time  until  the  next  occurrence.  Subdivide  the  time 
interval  of  any  fixed  length  t into  n equal  parts.  Given 
that  n is  sufficiently  large,  assumptions  ( 2 ) and  (3)  guarantee 
that  the  probability  of  at  least  one  occurrence  in  any  specific 
one  of  the  n parts  is  nearly  proportional  to  the  length  of 
the  interval  and  hence  can  be  written  as  Xt/n  for  some  constant  X>C. 
Then  the  probability  that  there  will  be  no  occurrence  in  any  specific 
one  of  the  n parts  is  1 -Xt/n.  The  random  variable  T is  greater 
than  t if  and  only  if  there  are  no  occurrences  in  any  of  the  n 
intervals.  Given  the  independence  affirmed  by  assumption  (1)  we  can 
then  write  the  approximate  equality, 

, t n 

P(T>  t)  * (1  - . 

n 

This  relationship  becomes  an  actual  equality  in  the  limit  as  the  length 
of  the  n parts  goes  to  zero.  Letting,  then,  n go  to  infinity 
we  obtain  by  a basic  limit  theorem, 

-Xt 

* e 


P(T  > t) 
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The  distribution  function  of  T is  then  given  by 

PT(t)  = P(  T < t)  = l-e~Xt  , t>0, 

and,  by  differentiation  with  respect  to  t, 

f (t)  - — FT(t)  = xe'U,  t > 0. 

1 at 

This  completes  the  derivation  of  the  asserted  form  of  the  proba- 
bility density  function  of  T. 

One  consequence  of  assumption  (3)  is  that  the  rate  X at 
which  dyadic  wars  are  being  initiated  does  not  change  over  time.  This 
could  pose  difficulties  in  the  present  research  in  two  respects:  yearly 
and  secular  variations.  In  the  first  respect,  Richardson  (i960,  p.  129) 
notes  a tendency  for  wars  to  begin  in  the  spring  and  fall,  but  since 
years  are  used  as  the  time  period,  for  the  aggregated  aata,  this 
effect  is  removed.  The  second  problem  of  secular  variations  across  a 
number  of  years  is  more  difficult  and  does  pose  certain  problems  which 
will  be  discussed  in  Section  5.2. 

Assumption  (2)  is  violated  when  simultaneous  occurrence  of 
the  initiation  of  two  or  more  wars  is  possible.  Thus  it  is  illadvised 
to  model  the  number  of  nations  initiating  wars  using  the  Poisson  process 
since  any  initiation  is  necessarily  dyadic  and  would  violate  the  assump- 
tion. The  use  of  dyads,  however,  is  not  in  obvious  violation  of 
this  assumption. 

i.U  _i. 

Let  us  denote  the  time  between  the  i and  (i+l)55  dyadic 
war  initiation  in  the  period  1816-1965  as  Ti  for  i-1,2, . . . ,208. 

Then  we  can  consider  the  ordered  (in  time)  set  (T-^Tg,  • • • »T20g)  as 
a time  series  of  length  208.  As  stated  above  the  assumptions  of  the 
Poisson  process  guarantee  that  these  random  variables  should  be  independent 
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and  exponentially  distributed.  A test  of  the  Poisson  assumptions  can 
be  made  by  determining  if  these  random  variables  are  independent. 

A standard  method  for  testing  independence  of  elements  of  a time  series 
is  through  sample  autocorrelation  coefficients.  The  k lag 
(Pearson)  sample  autocorrelation  coefficient  r^  is  defined  by 

rk  = S(Tt  - T)(Ti+lc  - T)/S(T1  - T)2, 

where  T = 2 T,/n  and  n is  the  length  of  the  series.  If  long(  short) 

i 1 

waiting  times  are  immediately  followed  by  long( short)  waiting  times, 
the  first  (i.e.,  k=l)  lag  autocorrelation  coefficient  should  be  positive. 
Similar  remarks  apply  for  k=2,3,...  • 

When  the  time  series  is  assumed  to  be  normally  distributed, 
Pearson  lag  autocorrelations  are  usually  employed  to  test  independence. 
Since  the  waiting  time  series  is  thought  to  be  more  nearly  exponentially 
distributed,  we  instead  use  Spearman ‘lag  autocorrelations.  As  in  the 
familiar  test  of  independence  using  Spearman' s rho,  ranks  of  the 
waiting  times  are  employed,  making  this  approach  ron-parametric  (Cox 
and  Lewis,  19 66,  p.  166  ).  For  the  waiting  time  series,  the  value  of 
Spearman's  first  lag  autocorrelation  coefficient  is  t^s.269,  for 
which  the  p-value  is  < .001.  This  suggests  positive  correlation 
between  successive  waiting  times.  Similar  positive  Spearman  auto- 
correlation coefficients  are  found  for  k*2,3,4.  This  shows  that  the 
process  of  dyadic  warfare  initiation  is  not  consistent  with  randomness 
and  also  suggests  the  direction  of  the  deviations  from  randomness. 

Although  it  appears  unlikely  that  Poisson  assumptions  are 
satisfied,  for  comparative  purposes  we  report  a test  of  fit  of  the 
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number  of  war  initiations  per  year.  Let  denote  the  number  of 

dyadic  war  initiations  observed  in  the  year  1815 + i for  1*1,2, . . . ,150. 
Then  under  the  Poisson  assumptions,  X^,  X2,  • • • > X^Q  be  a sample 

of  size  150  from  a Poisson  distribution  with  parameter  By 

estimating  \ an  estimated  theoretical  distribution  can  be  obtained. 
Table  1 reports  the  observed  distribution,  and  the  theoretical 
distribution  implied  by  the  Poisson  assumptions.  An  inspection  of  the 
differences  clearly  reveals  that  the  fit  is  poor.  The  chi-square  test 

p 

of  goodness  of  fit  (x  *169.29,  DF  = 4,  p<.001)  confirms  this.  It 
is  evident  that  unlike  the  instances  of  war  (as  defined  by  Richardson 
(i960)  and  also  by  Singer  and  Small  (1972))  the  process  of  dyadic 
warfare  Initiation  cannot  be  adequately  modeled  by  the  Poisson  process. 
The  next  sections  investigate  the  adequacy  of  more  complex  stochastic 
models  for  dyadic  warfare  initiation  including  formulations  of  hetero- 
geneity and  of  contagion. 

Table  1 about  here 


5.2.  Heterogeneity  over  time 

As  was  noted  in  Section  2 the  observed  cyclical  patterns  of 
warfare  may  suggest  wave-like  variation  in  the  Poisson  rate  parameter. 
Certainly  this  is  a possibility  suggested  by  the  graph  displayed  as 
Figure  1. 

Figure  1 goes  about  here 

We  Investigate  this  possibility  in  the  following  way:  In  agree- 

ment with  the  analysis  of  Denton  (1966),  and  Denton  and  Phillips  (1967) 
a Poisson  model  Is  postulated  with  the  probability  of  an  event  In  the 
interval  t to  t + dt  proportional  to  an  oscillating  function  of  t 
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with  increasing  amplitude.  Consistent  with  this  let 

X(  t)  = (c  + dt) (cos  yt  + 1)  » 2(c  + dt)cos2(yt/2) , 

c,  d>0,  the  one  being  added  to  insure  \(t)>0.  If  Xi  is  as 
in  Section  5*1>  then  it  is  known  (Parzen,  1962,  p.  125  ) that 

under  these  assumptions  that  X.^  is  a Poisson  random  variable  with 
rate  parameter 

-i+l 

\(u)du. 

J i 

Carrying  out  the  integration  we  find  that  the  rate  is 

r 2 

2 ( c + ud ) co s ( y u/2 ) au  = 

J i 

p 

c + di  + d/2  + d(  cos  y(  i +1)  - cos  yi  )/y  + 

( ( c + d(  i+l) ) sin  y ( i+l)  - ( c + di)  sin  y i)/y  . 

Table  1 reports  the  results  of  fitting  this  model  of  the  data 

on  dyadic  wars:  The  agreement  between  the  observed  and  theoretical 

2 o 13 

distributions  is  poor  (x  = 39*06;  DF = 4;  p=7xl0~°).  This  suggests 
that  the  distribution  of  dyadic  wars  in  time  is  not  adequately  modelled 
by  the  Poisson  process  with  cosine  time  dependency.  An  alternative  to 
this  model  would  be  \(t)  = ec+dt[cos  Xt  + 1]  , but  this  is  unlikely 
to  make  any  substantial  difference. 

We  therefore  examine  some  alternative  models  consistent  with 
the  conceptual  framework  of  Sections  2 and  3. 
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50  Heterogeneity  over  actors 

As  stated  by  Taibleson  (1974,  p*  878)  "there  has  been 
consensual  agreement  that  the  assumption  of  heterogeneity  over 
actors  will  lead  to  the  Greenwood-Yule  model".  This  model  specifies 
a particular  form  of  heterogeneity  in  that  a number  of  Poisson 
processes  are  operative  ana  the  rate  parameters  are  not  necessarily 

14 

equal  and,  in  fact,  vary  according  to  a gamma  distribution.  Thus, 
according  to  this  model  the  individual  nations  in  the  system  have 
different  rates  of  participation  in  warring  dyads  (i.e.,  heterogeneity), 
unlike  the  Poisson  process  where  the  rates  are  uniform. 

The  Greenwood-Yule  model  leads  to  the  number  of  events 
in  a specific  time  interval  given  by  the  negative  binomial 
distribution.  If  the  parameters  of  the  gamma  distribution  are 
a and  g , the  probability  pk  of  obtaining  k:  dyadic  ware 

initiations  in  any  time  period  of  unit  length  is 

eV(f3+l)k+a  k=0, 1,2,...  . 

For  a derivation  of  this  fact  see  Johnson  and  Kotz  (1969,  P*  25). 

Moment  estimates  for  0 and  a are 

9 = ( s2  - x)/x  and  a = x2/(s2-x). 

In  this  example  the  estimates  are  a=  .1255  and  9 = 11.04.  The 
expected  distribution  is  shown  in  Table  1.  While  the  fit  of  this 
distribution  is  much  better  than  the  Poisson,  the  differences  are  still 
substantial  enough  so  that  statistically  the  differences  are 
significant  (x2- 19*98,  DF  * 4,  p«5xlO”^). 
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Thus,  although  allowing  for  different  rates  of  participation 
among  the  warring  dyads  gives  a better  fit,  the  model  is  still 
deficient. 

5.4  Addiction 

A general  model  to  represent  the  situation  of  addiction1^ 
is  to  assume  that  each  of  the  dyads  follow  independent  Poisson 
processes  each  with  a possibly  time-varying  parameter.  If  there  is 
positive  (negative)  addiction,  the  rate  of  a particular  dyad  should 
go  up  (down)  each  time  that  it  initiates  hostilities.  If  the  rates 
are  assumed  to  vary  according  to  \ = a+bj  where  j is  the  number 
of  previous  times  the  dyad  has  fought  and  a and  b are  parameters, 
the  model  is  called  the  Eggenberger/Polya.  This  model  allows  for 
the  possibility  of  both  positive  and  negative  addiction — depending  on 
the  sign  of  the  parameter  b. 

Arbous  and  Kerrick  (1951,  p.  411)  showed  that  the  Eggenberger/ 
Polya  model  leads  to  the  negative  binomial  distribution  for  the 
number  of  dyadic  wars  in  any  time  period.  The  parameters  a and 
9 of  the  negative  binomial  distribution  are  related  to  a and 
b by 

a * a/b  and  0 = eb  - 1. 

The  fit  of  the  Eggenberger/Polya  is  then  exactly  the  same  as  that 
of  the  Greenwood- Yule,  which  from  the  discussion  in  Section  5-3  we 
know  is  poor. 

Although  the  Polya/Eggenberger  is  a natural  contagious  model 
for  certain  fields,  such  as  accident  statistics  (Arbous  and  Kerrlck, 
1951  > Bates  and  Neyman,  1952  ),  its  restrictive  assumptions  make 
it  doubtful  that  it  can  be  used  to  explain  dyadic  warfare.  It  seems 
reasonable  that  the  rate  of  the  various  dyads  should  also  be  a function 
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of  the  number  of  wars  initiated  recently  by  other  dyads.  To 
investigate  this  conjecture,  we  new  turn  to  infectious  contagion 
models. 

5.5  Infection 

In  contrast  to  addiction,  in  an  infectious  process  the  rates 
of  all  dyads  may  change  when  some  dyad  goes  to  war.  If  the  rates  are 
increased  (decreased),  at  least  temporarily,  we  say  that  the  process 
has  positive  (negative)  infection.  Most  and  Starr  (1976)  call  this 
positive  (negative)  spatial  diffusion.  If  there  is  positive  infection, 
long( short)  times  between  the  start  of  dyadic  wars  should  be  followed 
on  the  average  by  long  (short)  times. 

To  statistically  investigate  the  possibility  of  infection 
we  use  a contingency  table  approach  based  on  the  time  between  successive 
entry  of  dyads  into  war.  A two-way  contingency  table  is  constructed  by 
picking  numbers  a , a2,...,ar  with  0 = aQ<a1<. . . <ar,  and  then  classifying 
the  point  (Ti,Ti+1)  as  belonging  in  cell  (j,k)  if  aj._1<Ti<a_  and 
ak-1 < T^+1 < a^.  An  Infectious  process  suggests  a certain  pattern  to 
the  cell  entries  in  the.  contingency  table:  If  there  is  positive  infection, 

most  of  the  pairs  should  fall  on  (or  near)  the  main  aiagonal  ( i.e.,  upper 
left  to  lower  right).  If  there  is  negative  infection,  most  of  the  pairs 
should  fall  on  the  transverse  diagonal  ( i.e.,  upper  right  to  lower  left).1^ 
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Ths  Poly bs r model  Implies  that  the  dixierence 
In  time  between  successive  wars  should  be  approximately  independent  and 
should  have  an  exponential  distribution.  Thus,  each  of  the  rows  of  the 
contingency  table  should  have  approximately  the  same  frequency 
distribution. 

In  order  to  make  the  cell  entries  large  and  to  have  the 
same  number  in  each  row,  a 3x3  contingency  table  was  used  with 
divisions  0,  1-100,  and  >L00  ( in  days).  Table  2 shows  the  number 
observed  in  the  3 intervals,  the  number  expected  under  independence, 
and  the  difference  between  these  numbers. 

Table  2 about  here 
2 

The  standard  x test  indicates  deviations  from  independence 
(x  =11.3;  DF  = ^;  p=.02).  Furthermore,  the  differences  strongly 
suggest  the  deviations  are  in  the  direction  of  positive  infection.  The 
differences  on  the  diagonal  are  all  positive,  and  those  off  the  diagonal 
are  typically  negative.  Similar  conclusions  would  be  reached  with 
contingency  tables  of  other  dimensions  and  other  division  points. 
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There  is  no  reason  to  assume  that  only  the  length  of 
time  between  the  two  most  recent  wars  influences  the  Foisson  parameters. 
If  these  parameters  are  influenced  by  the  difference  in  time  between 

^ ^ 4- 

the  (j-l)3  and  j previous  outbreak  of  war,  there  should  be 
correlation  between  and  T^+j.  We  use  the  same  contingency  table 

analysis  to  test  for  independence  and,  more  importantly,  to  find  the 
direction  of  the  deviations  from  independence.  A chi-square  test  of 
independence  using  the  same  methodology  as  above  was  carried  out 
for  J»2,3»  • • • >9-  The  chi-square  and  p-values  are  shown  in  Table  5. 


Table  3 about  here 


These  tests  show  that  more  than  just  the  time  between 
the  last  two  wars  are  important.  The  chi-square  statistics  are 
significant  at  the  5^  level  for  j=l,2,4,5,  and  7.  Furthermore, 
the  same  pattern,  characteristic  of  positive  infection,  is  evident  in 
each  case.  The  dependence  dies  out  as  j increases. 

From  this  analysis  it  is  evident  that  some  positive  infection 
is  present  in  the  system  and  that  more  work  is  necessary  in  modelling 
it  in  order  to  more  precisely  capture  the  character  of  the  process. 

In  the  next  section  a parametric  model  is  developed  which  provides  a 
quite  acceptable  fit  to  the  data. 

5.6  An  Infectious  Model 

A model  that  allows  for  infection  is  the  autoregressive 
process  of  order  p (AR(p))  given  by 

\ “Jj  * i Ht-i+  5 + V 
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where  jz^, • • • »0pJ6 • and  P are  parameters  and  the  e^  are  independent 
and  identically  distributed  (i.i.d. ) random  variables  each  with 
mean  O.1^  We  model  the  interarrival  times  as  an  AR(p^  process  • 

after  using  the  logarithmic  transformation  = log(Tk+.5)*  The 
logarithmic  transformation  is  known  to  be  beneficial  in  modelling 
series  when  the  percentage  change  is  more  homogeneous  than  the 
actual  change.  ° If  the  k n waiting  time  Tk  was  recorded  as  0 
lays,  the  actual  time  was  some  fraction  of  a day.  Thus,  a better 
approximation  (on  the  average)  to  the  true  waiting  time  is  obtainea 
by  using  .5  instead  of  0.  Since  the  transformation  log(T+c)  is 
only  sensitive  to  the  value  of  c (0<c<l)  for  small  T,  we  use 
c=.5*  Thus,  the  assumed  model  is 

4 

log  (Tk+.5)  = 2^  ^ log(Tk_i+ .5) +6  + ek.  (2) 

In  the  case  p = 0,  the  model  (2)  can  be  written 

log(Tk+  .5)  = 6 +ak,  (5) 

where  a^  are  i.i.d.  random  variables  with  mean  0.  This  model 
implies  that  the  variables  log(Tlc+.5)>  hence  Tk,  are  i.i.d. 
random  variables.  Since  the  Poisson  assumptions  imply  that  the  waiting 
times  T^  are  i.i.d.,  the  AR(0)  process  agrees  with  the  Poisson 
assumptions.  Although  the  qualitative  statement  that  the  Poisson 
assumptions  do  not  hold  has  been  clearly  demonstrated,  introduction 
of  model  (J>)  is  useful  in  determining  the  amount  of  infection  in  the 
system  since  it  represents  a baseline  case. 

The  order  p of  the  autoregressive  process  can  be  estimated 
through  the  partial  autocorrelation  function.  If  the  model  is 
actually  AR(p),  for  J > p the  J th  order  (estimated)  partial 
autocorrelation  r\  < is  approximately  normally  distributed  with  mean 
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0 and  variance  n . Thus,  for  example,  for  j > p 

P(hj|  > 2.85n‘^)  = .004. 

Suppose  we  take  the  reasonable  position  that  we  are  willing  to  consider 

only  p such  that  p<m,  where  an  order  larger  than  m would  not 

be  practical  or  lack  plausibility.  The  backward  method  of  selecting 

the  order  of  the  process  is  accomplished  by  picking  positive  constants 

Cj  for  1<  J<m  and  then  picking  p to  be  the  largest  J such 

that  J ti  j | >Cj.  The  constants  Cj  are  picked  such  that  if  model  (3) 

is  true,  the  probability  of  obtaining  any  significant  partial  auto- 

m 

correlation  (l  - II  P(|r|.|  < c.,  ))  is  small.  Suppose  we  pick 

J J 

m=l2  and  c^  = 2.85n  ^ for  1<^^12.  Then  the  level  of 
significance  of  the  backward  method  is  approximately  a = 1 - ( -996)1*  = . 05- 
The  first  12  partial  autocorrelations  for  the  inter- 
arrival times  are  shown  in  Table  4.  Based  on  the  backward  method 

Table  4 about  here 


iiscussed  above,  we  estimate  p = 4.  The  least  squares  estimates  of 
the  parameters  of  model  (2)  '(using  natural  logarithms)  with  p = 4 are 

&L  = -269,  }S2  = .132,  ^ = .068,  2^  = .193,  and  6 = .777. 

The  fact  that  all  four  autoregressive  parameter  estimates  are  positive 
suggests  positive  infection.  That  is,  this  model  implies  long( short) 
interarrival  times  follow  long( short)  interarrival  times. 

The  above  analysis  suggests  that  there  is  some  advantage 
in  using  model  (2)  to  forecast  time  until  the  next  dyadic  war  initiation. 
To  quantify  the  additional  accuracy  obtained,  we  assume  that  the 
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forecasting  loss  function  is  given  by  square  error.  Table  6 gives 
the  minimum  mean  square  error  predictor  of  log(Tn+1+.5)  given  the 
previous  values  of  the  series  (i.e.,  for  l<i<n).  The  forecasts 
for  log(Tn+1  + . 5)  can  easily  be  converted  into  predictions  for  Tn+1 
by  exponentiation. 


residual 


The  forecast  accuracy  can  be  measured  through  the  model 

4 

s given  by  ek  = log(Tk+.5)  -o-Z  jzf^  log(Tk  i+-5)  for  model 


n 


(2)  and  ak  = log(Tk+.5)  -2  log(Ti+.5)/n  for  model  (3).  The  fore- 
casting error  of  the  next  value  of  log(Tn+1  + .5)  can  be  estimated 

n „ 

by  the  residual  mean  square  error  given  by  2 e£/(n-4)  for  model  2 

5 K 

n » 2 

and  by  2 a 7,/n  for  model  (3).  These  values  are  given  in  Table  5* 

1 K 

The  percentage  reduction  in  mean  square  error  obtained 
by  the  introduction  of  the  autoregressive  parameters  is  estimated 


by 


r2  = 1 - 5 


n 

2 e2/(n-4) 


n *2 
2 a5/n 
1 K 


Using  the  data  from  1816-1965  we  find  r2  = 1 - 7.105/9.163  = .225- 
Thus,  we  obtain  a 22.5  percent  reduction  in  forecasting  variance 
by  introducing  the  autoregressive  parameters.  This  22.5  percent 
reduction  gives  a numerical  measure  of  the  extent  to  which  infection 
is  present. 

A chi-square  goodness  of  fit  test  based  on  the  residual 

autocorrelation  function  is  available  (Box  and  Jenkins,  1970,  p.  291). 

If  rk  is  the  k lag  (Pearson)  sample  autocorrelation  coefficient  of 

the  residual  series  e.  (i.e.,  replacing  T^  with  e<  in  the  defining 

2 ® 2 

formula  for  rk  in  Section  5*1),  the  x statistic  is  Q ■ N I 


where  N is  the  number  of  residual  ej.  Using  m=30  the  value  of  Q 

2 

is  31.65  which  should  be  compared  with  a x2^  null  distribution. 
This  gives  a p-value  of  .17,  which  suggests  that  the  fit  of  the 
autoregressive  model  of  order  4 is 


reasonable. 
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6.  CONCLUSIONS  AND  IMPLICATIONS  FOR  FUTURE  RESEARCH 

The  basic  conclusion  of  this  article  is  that  the  process 
of  dyadic  warfare  initiation  is  predominantly  contagious,  and 
specifically  infectious.  This  is  a statement  descriptive  of  the 
stochastic  process  of  war  and  it  has  important  implications  for 
the  prediction  of  the  growth  of  war.  Specifically,  tne  model 
developed  in  Section  5*6  for  the  interarrival  times  can  be  used  to 
forecast  the  time  period  to  the  next  dyadic  war  initiation.  The  fit 
of  this  model  must  be  considered  surprisingly  good  in  view  of  the 
discussion  in  Section  2 of  the  potential  impact  of  war  on  inter- 
national system  dynamics. 

Such  empirical  predictions,  of  course,  should  be  interpreted 
with  considerable  caution  and  work  begun  on  a causal  interpretation 
of  warfare  in  which  relevant  causal  variables  would  be  identified 
and  quantified.  Consistent  with  this  and  our  finding  of  infec- 
tiousness we  have  a conceptual  framework  which  suggests  that  one 
question  for  further  research  is  why  some  wars  spread  and  others 
do  not.  Rather  than  asking  what  causes  wars,  we  are  asking  what 
makes  some  of  them  grow.1^  We  might  seek  causal  variables  which 
influence  this  growth  process  by  looking  to  previous  research  on  the 
contagion  of  social  behavior.  Two  variables  which  immediately  attract 
our  attention  are  interaction  opportunity  and  status.  The  basic  con- 
cept  here  is  one  of  centrality  (Reynolds,  1971).  Midlarsky  (1970) 
found  among  Latin  American  nations  a strong  relationship  between  a 
nation's  diplomatic  status  and  the  diffusion  of  military  coups.  A 
substantial  amount  of  research  in  other  areas  indicates  that  the 
behavior  of  high  status  individuals  is  associated  with  subsequent 
diffusion  of  that  behavior  (Rogers  and  Shoemaker,  1971). 
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There  is  also  research  which  supports  the  proposition 
tnat  contagion  within  a population  is  strongly  associated  with  the 
extent  of  interaction  within  the  population.  Coleman's  research  on 
the  diffusion,  or  contagion,  of  medical  innovations  indicates  quite 
clearly  that  physicians  who  had  the  greatest  interaction  with  each 
other  were  also  those  who  adopted  the  innovation  the  most  rapidly  and 
extensively  (Coleman,  Katz  and  Mengel  1957  ).  In  the  domain  of 
international  politics  there  are,  of  course,  numerous  sorts  of  interaction 
processes.  To  name  but  a few,  geographic  regions  (Russett  (1967)), 
foreign  trade  (Alker  and  Puchala  (1968)),  as  well  as  cultural  bonus 
and  alliances  (Deutsch  and  Singer,  1964)  all  furnish  interaction 
opportunities  and  patterns.  Alliances,  since  they  reflect  political 
commitments,  perhaps  provide  a template  for  the  future  growth  of  wars. 
Previous  research  has  focused  upon  the  relationship  between  various 
alliances  patterns  and  international  warfare  (Singer  and  Small,  1966 
and  1968).  This  research  nevertheless  has  not  directly  addressed 
the  probability  of  & war  spreading  within  an  alliance  network. 

Further  research  might  profitably  pursue  another  line 
of  inquiry  as  well.  The  data  used  in  the  present  study  represent  a 
rather  special  conceptualization  of  international  war.  It  contains 
a smaller  body  of  data  than  would  have  been  the  case  had  we 

chosen  to  use  the  complete  data  set  of  Singer  and  Small,  or  of  Wright 
or  Richardson.  The  use  of  these  data  sets  would  in  one  way  or  another 
led  us  astray  from  our  focus  upon  the  extent  to  which  violent  conflicts 
of  substantial  size  (at  least  1,000  fatalities)  occurred  In  the 
nation-state  system.  Richardson's  data  included  numerous 
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actors  such  as  ethnic,  cultural  or  political  groups  existing 

within  a nation,  and  consequently  the  wars  and  actors  listed  were 

not  exclusively  International.  Wright’s  list  also  had  the  characteristic 

of  Including  more  than  international  wars.  Other  listings  exist, 

but  these  also  tend  to  focus  upon  more  than  international  wars,  as 

well  as  being  limited  to  a more  contemporary  time-span. 

While  these  data  collections  were  of  limited  use  to  us  in 
the  present  study,  they  do  have  potential  for  extending  this  research. 

In  addition,  there  exist  other  lists  of  wars,  generally  dealing  with 
more  contemporary  phenomena,  which  could  be  useful  in  studying  the 
extent  to  which  local  conflicts  of  various  types  tend  to  spread  and 
at  what  rates.  For  example,  Kende's  (1971)  research  on  local  wars 
presents  not  only  a listing  but  a typology  as  well.  Such  a list  could 
be  helpful  in  identifying  the  processes  of  contagion. 


#?* 


Table  1 Instances  of  Dyadic  War  Initiation  and  Expected  Poirson, 
Poisson  with  Cosine,  and  Negative  Binomial  Distributions 


i number  of  instances  of  n,  number  of  Poisson  Poisson 
dyadic  war  initiation  years  Expected  with 

Cosine 


0 

95 

37.49 

Expected 

66.43 

1 

27 

51.99 

31.69 

2 

8 

36.05 

19.66 

3 

6 

16.66 

12.94 

4 

3 

5.78 

8.34 

. 5 

2 

1.63 

5-10 

> 6 

9 

.45 

5-90 

p 

chi-square  (x  ) 

I69.29 

39-06 

degrees  of  freedom 

4 

4 

x = 1.387  s2  = 16.70 


Negative 

Binomial 

Expected 

109.76 

12.63 

6.52 

4.24 

3.03 

2.30 

11.52 


19.98 

4 


Table  2.  Relationship  Between  Successive  Interarrival  Times 


Observed  Numbers  for  T^+1 


0 

1-100 

>100 

Observed 

0 

51 

21 

21 

Numbers 

1-100 

27 

18 

11 

for  Tj_ 

>100 

15 

17 

25 

Expected  Numbers  for  Ti+1  Unaer  Independence 
0 1-100 >100 


Expected  Numbers  0 

42.0 

25-3 

25.7 

for  Tj_  Under  1-100 

25-5 

15.2 

15.5 

Independence  >100 

25-7 

15-5 

15.8 

0 1-100 >100 


Difference  Between 

0 

9 

-4.3 

-4.7 

Observed  and 

1-100 

1.7 

2.8 

-4.5 

Expected 

>100 

-10.7 

1.5 

9-2 

1825  1845  1865  1885  1905  1925 
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Table  3.  Values  of  x2  tests  for  independence 

*123^56789 

X2  11.3  14.1  6.9  31.5  11.6  4.8  10.5  6.4  4.6 

p .02  .007  .14  2 x 10‘6  .02  .31  .03  *l8  .33 

Table  4.  Partial  Autocorrelation  Coefficients  for  log  (T^+.5) 

lei  2 3 4 5 6 7 8 9 10  11  12 

*1  .319  .167  .034  .290  .056  .063  .141  -.172  -.024  .106  .056  .024 

Table  5.  Optimal  Predictor  of  log( Tn+1  + . 5)  and  Mean  Square  Error 

^ Model  2 Model  3 

Optimal  Predictor  S+I  ^ log  (Tn+l-i  + ’ 2.059 


Forecasting  (mean  square) 

error  7-105  9-165 


NOTES 
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A review  of  Richardson' s efforts  in  this  area  will  not  be  pre- 
sented here  since  an  extensive  description  of  the  twelve  models  is 
presented  in  Zinnes  (1975  and  forthcoming).  A more  general  evaluation 

of  the  Statistics  of  Deadly  Quarrels  may  be  found  in  Wilkinson  (1976). 

2 

For  a speculative  discussion  of  this  see  Duncan  (1976). 

^ When  the  data  are  divided  into  Nineteenth  and  Twentieth  Century 
periods,  the  resulting  tables  still  do  not  suggest  dependence. 

4 

The  Poisson  process  and  its  implications  will  be  discussed  in 
Section  5*1* 

Although  the  Statistics  of  Deadly  Quarrels  was  published  in  i960, 
much  of  the  research  it  contains  was  completed  between  1941  and  19^9; 
Richardson  died  in  1953 • 

g 

Other  sources  of  time  heterogeneity,  which  arise  from  the  feed- 
backs between  the  phenomenon  under  study  and  other  social,  economic, 
military,  or  political  processes,  are  not  so  easily  handled.  Thus, 
for  example,  in  a study  of  a long  term  conflict  environment  such  as 
the  Mideast,  it  must  be  realized  that  a major  event  such  as  a war 
may  have  serious  and  profound  consequences  for  the  social,  economic, 
military,  and  political  processes  of  the  system.  The  fact  that  warfare 
changes  the  factors  which  determine  the  future  incidence  of  warfare 
means  not  only  that  the  stochastic  process  in  operation  is  not  homo- 
geneous in  time,  but  that  actual  realizations  of  the  process  determine 
its  future  development.  Such  heterogeneities  cannot  be  dealt  with 
simply  by  letting  the  parameters  be  functions  of  time  and  resorting  to 
more  computationally,  involved  estimation  procedures.  Usually  the  only 
recourse  for  the  analysis  of  processes  this  complex  is  to  large-scale 
simulation  methods  on  a computer.  These  are  relatively  easy  to  imple- 
ment in  specific  situations. 

j 

Most  and  Starr  (1976),  writing  in  the  context  of  a study  of  the 
spread  of  warfare,  have  also  noted  the  importance  of  drawing  a dis- 
tinction between  these  two  types  of  contagion.  They  have  labelled 
a process  displaying  addiction,  a reinforcement  process,  while  a 
process  displaying  infection  is  called  a spatial  diffusion  process. 

Q 

Examples  of  this  are  to  be  found  in  Germany' s attacks  on  Czech- 
oslovakia in  1939  and  Denmark  in  1940,  as  well  as  the  annexation  of 
Austria  in  193° • Further,  the  Soviet  invasion  of  Estonia,  Lithuania 
and  Latvia  in  1939  found  no  resistance.  More  recently  India's  annex- 
ation of  Goa  was  not  resisted  by  Portugal. 

^ For  a discussion  of  Singer  and  Small's  criteria  for  including  an 
international  war  in  their  listing,  see  pp.  17-39  of  their  report  (1972). 

To  illustrate  the  potential  of  this  problem,  had  we  accepted  Singer 
and  Small's  dyads  it  would  have  necessitated  the  recording  of  the 
United  States  and  Finland  as  a warring  dyad  in  1942  when  Finland 
attacked  the  Soviet  Union.  However,  in  World  War  II  the  United  States 


and  Finland  not  only  failed  to  declare  war  on  each  other,  but  main- 
tained diplomatic  relations  until  1944. 


Judgments  as  to  who  actually  fought  whom  were  substantially  aided 
by  Richardson's  (i960)  matrix  of  each  war.  Wright's  (1965)  listing 
was  also  used,  but  inaccuracies  were  found  in  several  dates.  For 
example,  Wright  (p.  1557)  lists  Bulgaria's  declaration  of  war  against 
the  United  States  in  World  War  II  as  June  22,  1941,  almost  six  months 
prior  to  Pearl  HarborJ  Bulgaria' s actual  declaration  of  war  against 
the  United  States  took  place  on  December  13,  1941.  Further,  Wright 
recorded  many  members  of  the  British  Empire  as  entering  World  War  II 
in  January,  1942  when  in  fact  these  acts  were  taken  in  December,  1941 
(Royal  Institute  for  International  Affairs,  1947). 


12 

For  other  research  in  international  politics  which  relies  upon 
the  Poisson  process,  see,  e.g. , Job  (1973),  Siverson  and  Duncan  (1976), 
and  McGowan  and  Rood  (1975) • 


^ If  Y has  a chi-square  distribution  with  DF=4, 
can  be  obtained  from  the  equation. 

P(Y>x2)  * (1  + X2/2)e"x  /2. 


the  exact  p -value 
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The  gamma  distribution  specifies  a positive  random  variable,  say 

X,  whose  probability  density  function  is  proportional  to  xa-1e_x^ 
where  a and  0 are  positive  parameters. 

^ Most  and  Starr  call  this  phenomenon  reinforcement. 

An  analogous  procedure  to  one  suggested  by  Most  and  Starr  (1976) 
could  also  be  used.  That  is,  break  the  time  period  up  into  two  non- 
overlapping periods  and  for  each  dyad  record  the  number  of  times 
they  fought  in  the  two  periods.  These  observations  can  be  labelled 
(Ui,  V^),  i=l,...,n,  where  n Is  the  number  of  dyads.  A two-way 
contingency  table  analysis  could  be  used  on  the  pairs  (U.,V, ) as  in 
Most  and  Starr.  A problem  with  this  approach  is  that  many  Nations  do 
not  maintain  their  national  status  over  the  entire  time  record  because 
of  changes  in  national  boundaries,  annexations,  mergers,  etc.  This 
problem  becomes  more  severe  as  the  time  span  studied  increases  in 
length. 

17 

For  a more  complete  explanation  of  the  methodology  used  in  this 
section  see  Box  and  Jenkins  (1970)  or  Nelson  (1973)- 

18 

Acton  (1959;  P*  223)  has  written  "Data  that  are  counts  of  pop- 
ulations, vital  statistics,  census  data,  and  the  like  are  almost  always 
improved  by  taking  logs. .. Charles  Winsor  frequently  prescribed  the 
taking  of  logs  of  all  naturally  occurring  counts  (plus  one,  to  handle 
that  embarassing  quantity  zero)  before  analyzing  them--no  matter  what 
the  sources  [of  the  data]""!  Tufte  (1974;  p.  103)  has  suggested  several 
reasons  why  one  might  use  a logarithmic  transformation  in  regression 
studies;  an  Important  one  for  our  data  is  "Badly  skewed  distributions — 
in  which  many  of  the  observations  are  clustered  together  combined  with 
a few  outlying  values  on  the  scale  of  measurement — are  transformed  by 
taking  the  logarithm  of  the  measurement  so  that  the  clustered  values 
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are  spread  out  and  the  large  values  are  pulled  in  more  toward  the 
middle  of  the  distribution". 

This  it  may  be  argued  is  really  the  implicit  question  of  much 
research  on  the  causes  of  war.  Many  of  tne  variables  indicating  war 
are  closely  associated  with  the  size  of  war.  For  example.  Singer  and 
Small  (i960)  operationalize  war  as  a dependent  variable  by  measuring 
the  nation-months,  battle  casualties  and  the  number  of  wars  begun. 

The  first  two  are  reflective  of  size. 

20  Singer  and  Small  (1972)  have  presented  deta  indicating  that 
centrality  is  not  associated  with  war  participation.  Their  data, 
however,  are  of  a geographic  nature. 
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